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Abstract. Soft hadron production is described as a two-step process, where the interaction of the partonic 
constituents of the colliding hadrons leads to the production of intermediate subsystems (fireballs), which 
decay subsequently into hadrons. The weights of the various final states are derived from the corresponding 
phase-space factors modified by empirical transition elements. The results compare well with data at 
energies between particle production thresholds and ISR energies. Special emphasis is put on correlation 
data, which offer the opportunity to shed some light on the question whether particle production proceeds 
via fireballs or strings. 



1 Introduction 

QCD is assumed to be the theory of the strong inter- 
action. Soft hadron production, however, is a nonpertur- 
bative process which at present cannot be calculated by 
QCD. Thus, the understanding of soft hadron produc- 
tion is still based on phenomenological approaches like, 
e.g., the dual parton (DPM) |], the VENUS |] or the 
PYTHIA-LUND models. Within these approaches 

hadrons are considered as composite objects consisting of 
partons (quarks and gluons) the interaction of which is as- 
sumed to proceed in two steps. In a first step excited sub- 
systems, usually called strings, are produced which decay 
afterwards into stable particles and resonances. 

Ingredients of these models are structure functions, 
parton-parton cross sections and fragmentation functions. 
The structure functions of the interacting hadrons con- 
tain the information about the momentum distribution of 
partons. From the cross sections for parton-parton inter- 
actions the number of strings produced in a single hadron- 
hadron scattering can be deduced, and fragmentation func- 
tions describe the decay of strings into hadrons. These ap- 
proaches make maximal use of the information available 
from lepton-hadron and lepton-lepton interactions as well 
as from general properties of the scattering amplitude like 
unitarity and analyticity. 

On the other hand the first attempts to understand 
multiple hadron production in hadronic interactions were 
based on statistical considerations and the observation of 
excited intermediate subsystems called fireballs (FB) (see 
the reviews [0,|[^). A modern version of a thermodynam- 
ical approach can be found, e.g., in ref. pO[ |. 
The Rossendorf colhsion (ROC) model 



|15|,|16,17| is basically a statistical approach in deriving 
the relative contributions of the various final channels by 
calculating their statistical weights from the phase-space 
factors, which are, however, strongly modified by empiri- 



cal transition matrix elements. It aims at describing soft 
hadron production in the energy region between particle 
production thresholds and ISR energies. For this purpose 
the basic ingredients of hadron production models are re- 
formulated in such a way that they are applicable at low 
energies as well. The ROC model is built as a minimal 
approach in the sense that the number of parameters is 
restricted to the minimum necessary to well reproduce the 
main features of the available data. Special emphasis is 
put on the consideration of short-range correlations, be- 
cause the FB concept yields a natural explanation of the 
observed phenomena. 

The present paper is organized as follows. In sect, ^the 
basic features of the ROC model are explained. Section |^ 
contains a comparison of selected experimental data with 
model calculations. In order to show differences between 
string and FB approaches with regard to correlation data 
the ROC model is contrasted with the PYTHIA-LUND 
model 1^. By means of a special version of the ROC model, 
where FBs are degenerated into single hadrons, the im- 
portance of FBs for the description of correlation data is 
demonstrated. Conclusions are summarized in sect. ^ 



2 The model 

The basic idea of a statistical approach consists in the as- 
sumption that the probabilities of formation of the various 
final states are proportional to their statistical weights. 
This idea was implemented by Fermi [l^ fifty years ago, 
but his model turned out to be applicable only at rela- 
tively low energies. At higher energies it is no longer a 
good approximation to assume that the whole initial en- 
ergy is randomly distributed among the final particles. 
Particles with high transverse momentum, e.g., are pro- 
duced with extremely low probability indicating that the 
final states are dynamically linked with the initial state. 



2 



H. Miiller: Soft hadron production in pp interactions up to ISR energies 



In contrast to this early attempt the ROC model is 
based on the following modified statistical approach. In- 
stead of calculating the statistical weights from the whole 
phase-space the dynamics of the interaction is implemented 
in form of empirical functions which either suppress cer- 
tain regions of the phase-space or impose additional non- 
statistical weights. We define a channel a by the number n, 
masses nii and quantum numbers of the final particles. 
The relative probability of populating a channel a is calcu- 
lated as the product of the Lorentz-invariant phase-space 
factor dL„(s;a) with the square of an empirical matrix 
element A^, which describes the dynamics of the inter- 
action process. Here, s = denotes the square of the 
total energy with p being the total four-momentum. The 
phase-space factor is defined as the integral over the mo- 
menta of the final particles with energy and momentum 
conservation taken into account, 



dL„(s; a) = dL„(s; mi, ... , m„) 



2e» 



(1) 

where the four-momentum of the i-th particle is denoted 
by Pi = (ei,Pi) with p| = mf. For numerical calculations 



the S function in eq. (|^) is removed by introducing a new 
set of 3n — 4 variables to replace the 3n three-momentum 
components. It is reasonable to choose a set of variables, 
which reflects the underlying physical picture of the inter- 
action process. In fig. the reaction 



a + b^ FBi 



FB 



N 



hi 



hn 



between hadrons a and h resulting in the production of n 
particles is schematically depicted. First, 2 < TV < n in- 
termediate particle groups called FBs are produced, which 
decay into so-called primary particles. The primary par- 
ticles define the channels for which the weights (||) are 
calculated. Among them are resonances, which decay sub- 
sequently into stable hadrons. The phase-space factor cor- 
responding to the diagram of fig. |] can be calculated ac- 
cording to ref. |pj| 



dL„ (s; oln) 



N 



JldMjLnAMi;ai) 



.1=1 



dLN{s;Ah,...,MN) (2) 
with the invariant masses of the FBs equal to Mj = 



\fPj = Y (Srii^'j) ^be final channel defined as the 
vector OLtq = (ai, . . . , ajv) of the decay channels of the in- 
dividual FBs. The probability of populating the channel 
a AT is given by 

dVF(s; a^v) oc dL„(s; ckat) . (3) 

Here, the square of the matrix element contains the 
dynamical input and is split into factors 

which describe the interaction A? resulting in the produc- 
tion of iV FBs, the production of hadrons A^^^ via the cre- 
ation of quark-anti-quark {qq) pairs , the invariant-mass 
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Fig. 1. Phase-space decomposition of a two-step process. A'^ 
FBs with masses Mi{I = 1 . . . TV) are produced in the inter- 
action of two hadrons a and b with four-momenta pa and pi,. 
The FBs decay in the second step, where the J-th FB disin- 
tegrates into ni primary particles with X//!=i ~ Reso- 
nances among the primary particles decay afterwards into sta- 
ble hadrons. A possible decay chain is shown. 



distribution of the FBs A^^ , the transverse A^ and longitu- 
dinal momentum distribution of the FBs, and, finally, 
some additional factors A^^ necessary for the calculation of 
the statistical weights. In the following subsections these 
factors will be discussed in more detail. 



2.1 The interaction A? 

In models like DPM § or VENUS ||l the colliding hadrons 
are considered as extended and composite objects con- 
sisting of an indefinite number of partons, the interac- 
tion of which is assumed to proceed via color exchange. 
If, 6.17., the color exchange between the valence quarks of 
two baryons takes place, then two strings are produced. 
Each of them consists of the remaining diquark and the 
valence quark removed from the other baryon. Since more 
complicated exchanges are possible and the collision may 
proceed at different impact parameters a varying number 
of strings is produced. The corresponding probabilities are 
derived from Gribov-Regge theory ||2^ in the limit of high 
energies in combination with the use of profile functions 
for integrating over the impact parameter. PYTHIA |^ 
describes low-pt events on the basis of the multiple inter- 
action model of ref. Q , which extends a high-pt picture 
down into the low-pt region by regularizing the pt scale. 
The number of (independent) parton-parton collisions in 
one event depends on the impact parameter and on the as- 
sumed matter distribution inside the interacting hadrons. 

In the ROC model, non-statistical weights A? (TV) for 
producing a definite number TV of FBs are introduced. 
This is a phenomcnological parameterization of the con- 
tributions from the different color exchange diagrams and 
of the integral over the impact parameter, which we apply 
also at low energies. Such a dynamical input is necessary, 
because the phase-space factor alone tends to overestimate 
the number of FBs simply due to the fact that the number 
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of states increases tlie more FBs are produced. A thermo- 
dynamical approach using a chemical potential as the only 
parameter to regulate the number of FBs turned out to 
be not flexible enough. Therefore, the negative binomial 
distribution is applied where we have two free parameters 
V and q 

AnN;q,v)^(^~^y-qfil-qr. (4) 

In the calculations we use the mean N and the ratio of 
the variance to the mean D = a'^ /N as parameters from 
which q^[D~l)/D and v = N/{D - 1) follow. By this 
means, the whole complicated interaction scenario is de- 
scribed by altogether two parameters in conjunction with 
the corresponding phase-space factors. 

2.2 Quark statistics Al^ 

The factor stands symbolically for the algorithm ap- 
plied to sample the possible final states. At first the va- 
lence quarks of the interacting hadrons are redistributed 
among the FBs and then the final hadrons are produced 
via the creation of new qq pairs. All internal quantum 
numbers are conserved automatically by this procedure. 

The complicated details of color exchange diagrams 
and string drawings are replaced by statistical consider- 
ations. For describing the interaction process the notion 
of quark removal is borrowed from ref. It is assumed 
that the multiple interaction of the partonic constituents 
of the incoming hadrons leads to the creation of N FBs. 
Color exchange between the constituents results in the 
removal of the involved quarks and gluons from the in- 
coming hadrons. Removed quarks are found in any of the 
other FBs with equal probability. The remaining partons 
of the interacting hadrons form the two leading FBs, the 
scattered partons the N — 2 central FBs. In a next step 
qq pairs are produced and randomly (the q's and q's inde- 
pendently) distributed between the FBs such that each FB 
becomes color neutral and contains the minimal number of 
g's and g's necessary for building at least one hadron (me- 
son or baryon). This procedure is the equivalent for the 
sum over the possible color exchange diagrams with the 
restriction that the removal of two or three valence quarks 
or of a single sea quark is neglected. Only the removal of 
one valence quark with probability Wv is considered. The 
remaining probability 1 — Wv is understood as gluon or 
quark-pair removal. As equivalent to string fragmentation 
each FB is then filled separately with an arbitrary num- 
ber Uqi of additional qq pairs where nqj > 0. Up, down, 
strange and charm quarks are produced in the ratios 

M : d : s : c = 1 : 1 : As : Ac (5) 

with As and Ac being suppression factors due to the heavier 
masses of the strange and charm quarks. The creation 
of top quarks can be neglected in the considered energy 
range. 

The final hadrons are built up in each FB indepen- 
dently according to the rules of quark statistics |£lj by 



randomly selecting sequences of g's and g's. A qq gives 
a meson, while baryons or antibaryons are formed from 
qqq or qqq. From a given sequence of quarks the differ- 
ent hadrons are formed according to the tables of the 
particle data group ||2^. All baryons marked in the ta- 
bles with three or four stars, the meson nonets built from 
u, d and s quarks with angular momenta zero and one 
(^ 5*0, 5*1,^ -Pi,^ Pq,^ Pi,^ P2) as well as all charmed hadrons 
are taken into account. An empirical probability distribu- 
tion 

Wh(mh) cx exp (-TOh/Qh) (6) 

with an adjustable parameter 0h is used to suppress the 
formation of the heavier hadrons of mass rrih- During event 
generation the current masses of resonances are sampled 
according to a probability distribution consisting of the 
product of a relativistic Breit-Wigner distribution, the 
phase-space factor of the decay products and the above 
suppression factor Wh ^mh). The decay of resonances into 
the various channels proceeds either according to known 
probabilities or in accordance with the statistical weights 
of the possible final states in case of unknown decay prob- 
abilities. 

The described algorithm together with the parame- 
ters As, Ac and (9h is the equivalent of the usually much 
larger number of parameters describing the fragmentation 
of strings. 

It should be stressed that the ROC model has no pa- 
rameter fixing the probability of diffractive processes. A 
diffractive process is usually assumed to proceed via the 
exchange of a Pomeron, a fictitious particle which does 
not affect the quantum numbers of the involved particles. 
In the string models diffractive scattering is treated 
as a special process whose probability is determined by a 
free parameter adapted to data. In the ROC model diffrac- 
tive scattering is one of the possible final channels, because 
there is a certain probability that one or even both leading 
FBs are identical with the initial protons. This happens 
if the valence quark content of the considered FB remains 
unchanged, if no additional qq pairs are produced, and if 
a proton, and not a resonance, is built from the available 
quarks uud in the recombination phase. 

2.3 Mass distribution of FBs A^^ 

Until now we have explained how the hadrons forming the 
final state are sampled. In a next step the integral over the 
invariant masses of the FBs Mi [see eq. (||)] is performed. 
Again phase-space alone produces too large FB masses 
because of their corresponding large numbers of states. 
To restrict the invariant masses the FBs are assumed to 
be characterized by a temperature O. As matrix element 
squared the function 

N 

Al,^Yl{Mi/0)K^{MilO) (7) 

7=1 

with the asymptotic behavior 
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{Mi/0)Ki{Mi/O) 



1 



2.5 Statistics A^^. 



is used. The expression (Mr /0)K^ (Mr 10) is the kernel 
of the so-cahed _ft'-transformation (see used to trans- 
form a micro-canonical phase-space distribution depend- 
ing on the total energy Mi of the /-th FB into a canonical 
one, which is characterized by a temperature 0. In eq. 
Ki stands for the modified Bessel function. For increas- 
ing invariant masses M/ the function {Mj /0)Ki{Mj /0) 
strongly decreases, while the phase-space factor Ln,{Mi) 
of the /-th FB becomes larger. Their product has a max- 
imum at a value of Mj determined by the parameter 0, 
which thus fixes the average internal excitation energy of 
the /-th FB. Since FBs consist of only a few particles, we 
are far from the thermodynamical limit. Therefore, the 
momenta of the hadrons are calculated from the decay of 
the FBs according to phase-space and not from macro- 
canonical distributions. 



2.4 Longitudinal phase-space of FBs and 

FBs are proposed to emerge from the interaction of par- 
tons, whose momenta inside fast moving hadrons are mainly 
longitudinal. Due to the finite size of hadrons and the un- 
certainty principle a small transverse component is present 
too. This is taken into account by damping large trans- 
verse momenta of the FBs using a linear exponential dis- 
tribution 



In conclusion, all factors still necessary for a correct count- 
ing of the final states are collected in the term 



N 



= H exp(-7Pt,i) 



(8) 



1=1 



with the mean Pt = 2/7. The parameter 7 used here can 
be directly compared with the analogous parameter em- 
ployed in string models There, the transverse mo- 
menta of the partons, which form the ends of the strings, 
are restricted, here, we constrain the transverse momenta 
of the FBs directly. 

The longitudinal momentum distribution of the two 
leading FBs is weighted by 



Af 



with the scaling variables 



(9) 



Xi = {El + Pz,i)/\/i and X2 = {E2 - Pz,2)/Vi . 

Here, it is assumed that FB 1 is the remnant of the in- 
coming hadron a moving in the positive z-direction prior 
to the interaction, while FB 2 stems from hadron h mov- 
ing in the opposite direction. On the average the leading 
FBs carry the largest part of the longitudinal momenta as 
a consequence of the weighting (^. This forces the other 
FBs to have accordingly less longitudinal momenta. In this 
manner the factor A^ is the equivalent of the structure 
functions used in refs. 

The method of calculatin g th e longitudinal phase-space 
of N FBs is taken from ref. |2^ with appropriate modifi- 
cations due to the presence of A^ . 



71 J —1 




N-l 



(10) 



It contains the spin degeneracy factors (2(7^ + 1), the vol- 
ume V in which the particles are produced with V = 
AttR^ /3 determined by the radius parameter R. The quan- 
tity g{ai) — (jjf^n/j^^ is the degeneracy factor for 

groups of Ufj identical particles in the final state of the 
/-th FB and prevents multiple counting of identical states. 



2.6 DifFerentia! cross section 

Summarizing the above considerations we define the num- 
ber of states in the decay channel aj of one FB 



d2i{ai) ^ g{ai) 



{2n) 



Y[i2a, + l)2m, 



dMi{^)Ki(^]dL„AMi;ai) (11) 



and the analogous number for the set of FB states 

N-l ( N 



dZAf(s) 



V 



(27r)3 



n2A//exp(-7Pt,i)Af(iV) 



.1=1 



{XiX2fdLMis;A'h,...,MN). 



(12) 



From eqs. (^Tj) and (|T^) a compact expression for the prob- 
ability of populating the channel otN = (ai, . . . , cin) [see 
eq. (H)] can be derived 



dW{.s;aN) = |ndZ,(ai)|dZjv(s) 

The corresponding cross section is written as 

dW{s; ckat) 



(13) 



da{s) = crin(s) 



Ejv Ea„ / dVK(s; a^v) 



(14) 



where the inelastic cross section crin(s) of the considered 
reaction serves as normalization. Any physical quantity 
of interest can be derived from eq. dl4) by summing the 
contributions from all channels and integrating over the 
unobserved variables 
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Fig. 2. The calculated FB distribution da/dN, eq. (dots 
in a) and the measured [ 2^ and calculated multiplicity distri- 
bution du/dn of charged particles n (dots and histogram re- 
spectively in b) from pp interactions at ^/s = 63 GeV. In a 
the dependence of the matrix element Af{N) (line) and of the 
FB distribution without A?{N) (histogram) are depicted sep- 
arately (these two curves are arbitrarily normalized). 



2.7 Adjusting parameters 

The parameters introduced in the previous subsections de- 
termine definite features of the production process. So the 
radius parameter is responsible for the multiplicity n of fi- 
nal particles. In eq. (^0|) the factor V^~^ appears implying 
^3(n-i)^ and increasing R in turn means that states with 
large number n of produced particles get higher weights. 
The temperature and the radius R determine the mean 
invariant mass of the FBs, via the excitation energy, R 
via the number of particles in the FB. Abundances of dif- 
ferent hadron species are fixed by the probabilities As and 
Ac for the creation of the different quark flavors. The pa- 
rameter ©h in eq. (|^) influences the selection of hadrons 
consisting of the same valence quarks. Finally, the mo- 
menta of the primary hadrons are composed from the su- 
perposition of the FB momenta, determined by /3 and 7, 
with the relative velocities of hadrons in the rest frames 
of the FBs, influenced by 0. There are correlations be- 
tween the parameters, and several parameter sets have 
been found giving similar results. 

In searching for a suited set of parameters we start 
at the highest energy ^/s = 63 GeV with a guess for all 
parameters except N and D, which are responsible for the 
number of produced FBs [see eq.(||)]. These parameters 
are varied until the multiplicity distribution of charged 
particles is reasonably reproduced. 

In fig. H the influence of the factor A?{N) on the FB dis- 
tribution 



dais) EcyJdWis;^^) 

^ O-in(s); 



diV 



N 



(15) 



is demonstrated. Without the factor Af{N) the cross sec- 
tion for producing FBs (the histogram in fig. ^) reaches 



its maximum at values much too high for reproducing the 
multiplicity distribution of charged particles. By combin- 
ing the increasing phase-space factor with the decreas- 
ing function A? (N) (the line in fig. ||a) a FB distribution 
is produced which reproduces the distribution of charged 
particles (fig. |^b) . The similarity of FB and particle mul- 
tiplicity distribution is due to the fact that on the average 
each FB emits the same number of particles. 

In order to describe data also at lower energies the pa- 
rameters Wv, 7, P and R were made energy dependent, 
while the other ones are kept constant. For the prob- 
ability Wv of valence quark removal the energy depen- 
dence follows from the consideration of the binary reac- 
tions pp nA'^'^ and pp Z\°Zi++ for which data in a 
wide energy region are available (see fig. In the present 
approach contributions to these reactions come alone from 
events with N = 2 FBs. A removed quark from one pro- 
ton belongs, after the interaction, to the other one and 
vice versa. For the reaction to proceed an exchange of an 
u with a d quark must happen according to 

p + p ^ uud + uud udd + uuu n + A'^^ . 

The calculated cross sections become zero for Wv = 0, 
while for = 1 the data are well reproduced at low, but 
overestimated by orders of magnitude at high energies (see 
the dashed curves in fig. p|). Therefore, the function 



1-W, = Ei,/{Wi + Sk) 



(16) 



with the excess energy — ^~2mp and one parameter 
Wx is used to describe the energy dependence. It provides 
a smooth transition between the two extremes, Wv = 1 
for excess energy — Q and Wv — > for large excess 
energies. In this way the data are reproduced quite well 
with the correspondingly adjusted parameter Wi (see ta- 
ble |l|) . Such an energy dependence is in agreement with the 



10 - 




Fig. 3. Cross section of the reactions pp nA^^ and 
pp — » A*^ as a function of the laboratory momentum. 
Experimental data (dots) from are compared with ROC 
model results (solid lines). The dashed curves are calculated 
with the probability for valence quark removal Wv = 1. 
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Table 1. The parameter set used for the calculations. ME 
stands for the matrix element squared. 



10' 



ME 


Parameter 


Parameter 


Equation 


'data 


Af 


iV=1.8 


D=1.4 














As=0.15 


Ac=0.05 












6>h = 250 MeV 


Wi=3 GeV 






16 
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= 300 MeV 


R=1.2...1.29 fm 




I 


ic 




1 10 


^? 


70 = 3.8(GeV/c)-i 


71 = 1 GeV 




I 






E 


A? 


/3o = 2 


/3i = 1 GeV 




9 


18 







parton picture of hadrons. Roughly speaking the number 
of partons increases with energy and, consequently, the 
probability for a valence quark being involved in the color 
exchange diminishes. 

The energy dependence of the parameters 7 and /3 [see 
eqs. (|^), (^)] cannot be determined in such a clean way. 
Still the best indication comes from the proton spectra 
at different energies. In a diffractive event at least one of 
the leading FBs is identical with the corresponding initial 
hadron. The additional folding with the momentum dis- 
tribution from the FB decay is absent and the transverse 
momentum dep ende nce of the diffractive peak (see figs. ^ 
and ^ in sect. |3.2| ) is directly governed by 7 and (3. A 
reasonable description at the different energies is achieved 
using 

7 = 70^5/(71 + ^5) (17) 



and 



(18) 



The parameter R is then adjusted to reproduce the 
mean multiplicity of charged particles (rich) at the differ- 
ent energies. Below incident momenta of lOOGeV/c the 
constant value of i? = 1.29 fm is used, because (rich) ceases 
to be sensitive to the value of R. With increasing ener- 
gies the value of R decreases smoothly towards 1.2 fm. All 
ROC model results in this paper are calculated with the 
parameter set summarized in table |^. 



3 Comparison with data 

The ROC model is implemented as a Monte-Carlo gen- 
erator which samples complete events from which nearly 
all kinds of measurable quantities can be deduced and 
compared with experimental results. The overall agree- 
ment between data and ROC calculations is quite good in 
the whole considered energy range between about ^/s « 
2.2 GeV and the highest ISR energy of = 63 GeV. In 
the following a few selected data sets concerning hadron 
multiplicities and the dependencies of differential cross 
sections on longitudinal and transverse variables are pre- 
sented. Then correlations are discussed more thoroughly, 
because they are sensitive to differences between string 
and FB models. 



10 ■= 



p + p Vs = 27.5 GeV 
D° A- 1- r r*K° r U K"°K"'K'" K° p- n 



D^D-D^AA*** K'- A p A''A**K" r? p" w p 
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10 1 10 

multiplicity (model) 



Fig. 4. Hadron multiplicities for pp collisions at = 
27.5 GeV. The experimental average multiplicities p5|, p6| , |27| 
are plotted versus the calculated ones. The dashed line indi- 
cates coincidence between data and calculations. Well repro- 
duced data tend to lie near this line. A* stands for yl(1520). 



3.1 Multiplicities 

Hadron abundances are determined by the strange and 
charm suppression factors As and Ac, respectively, and the 
hadron temperature 0h- There is no parameter which gov- 
erns the production of baryons directly. Instead, the algo- 
rithm for building hadrons from quarks is responsible for 
baryon production. The more quarks are available in the 
FB considered, the more probable it is to select at random 
a sequence of three q^s (or ^'s) from which a baryon (or 
antibaryon) can be formed. Since the FBs become bigger 
with increasing temperature and radius i?, these are 
the parameters, which determine how many baryons are 
produced. A rather complete data set of hadron multi- 
plicities for pp collisions at ^/s = 27.5 GeV is compared 
with model calculations in fig. ^. The agreement is quite 
impressive except for the yl(1520) where the deviation is 
rather large. An obvious reason for this discrepancy has 
not been found. 

Figure ^ should be compared with the result from the 
thermodynamical model of Becattini (see fig. 5 in ref. JTo[). 
There, three parameters, temperature, volume and sup- 
pression factor As are fitted, which have similar meanings 
as in our study, although the parameter values differ con- 
siderably from those used here. Nevertheless, the achieved 
accuracy of the description is comparable. We have, con- 
trary to ref. [|lO| , charmed particles included in fig. ^, and 
the suppression factor Ac is adjusted to reproduce these 
data points. 

In fig. ^ the energy dependence of mean multiplici- 
ties of different particle species is plotted. The mean mul- 
tiplicity of charged particles (rich) is mainly affected by 
the temperature O and the radius R. In order to repro- 
duce (rich) the radius R is adjusted at each energy. The 
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p + p hadron + X 




10 

Vs (GeV) 

Fig. 5. The energy dependence of mean multiplicities of 
charged particles (ch) , -k^ , tt" , , K~ , p and p for pp scatter- 
ing. Full and dotted lines are fits to the data |28[, the symbols 
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Fig. 6. The topological cross sections cr„ for the production of 
n charged particles as a function of the laboratory momentum. 
Dots are data [E5|, full lines ROC results. 



adapted values of R change by less than 10% in the con- 
sidered energy region. Since charged particles are mainly 
pions the energy dependences of the number of charged 
particles and of pions are similar. The slightly larger val- 
ues of TT"*" compared to 7r~ arise from charge conservation. 
Strangeness suppression causes a large gap between pi- 
ons and kaons. It should be noted that, e.g., the VENUS 
model reproduces this gap (see fig. 10.2 of ref. |^) with 
a suppression u : d : s = 0.43 : 0.43 : 0.14 differing by 
a factor of nearly two from the value used here, namely 
u : (i : s = 1 : 1 : 0.15. A possible explanation for this 
difference might be the mass factor in eq. (|l^). Strange 
particles have heavier masses what gives them a larger rel- 
ative weight. On the other hand, the phase-space factor 
becomes smaller the heavier the produced particles are. 
Hence, a direct comparison of the suppression factors as 
used in string fragmentation and in the present approach 
seems to be difficult. The production of antiprotons is sup- 
pressed by about two orders of magnitude, because only 
FBs consisting of at least three quark pairs have a certain 
probability to create a baryon-antibaryon pair. The direct 
influence of the phase-space factor becomes especially im- 
portant in the threshold region where the multiplicities are 
small. The different curves for pions, kaons and antipro- 
tons are well reproduced due to the phase-space factor, 
which strongly decreases with diminishing kinetic energy 
still available after particle production. 

After having reproduced the mean multiplicities of char; 
particles by adapting the radius parameter we consider the 
energy dependence of topological cross sections in fig. || 
without any further parameter adjustment. The typical 
behavior of the cross sections with a weak maximum and 
subsequent slow decrease for low and a continuous increase 
in the considered energy region for high multiplicities is 
well reproduced by our model. 
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are data |^^, histograms represent ROC results. Data and cal- 
culations are corrected for acceptance. For the data n is the 
observed multiplicity, while for the calculations n is the (true) 
multiplicity in the pseudo-rapidity range 77 < 4. 
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3.2 Differential cross sections 

In a next step the distribution of particles in phase-space 
is considered. We start with the charged particle density 
for various intervals of the observed multiplicity shown in 
fig. 0. Charged particles were measured in the pseudo- 
rapidity region |?7| < 4. In the calculations the geometrical 
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Fig. 9. Invariant differential cross sections as a function of 
Feynman's variable xf at various values of p^ as indicated in 
the figure. Symbols are data , histograms ROC results. The 
spectra on the right hand side are multiplied by factors of 5 
and 3, respectively. 



acceptance of the apparatus has been taken into account 
according to the curve given in fig. 2 of ref. ||2^ . The differ- 
ences in the multiplicities observed in the experiment and 
the multiplicities in the calculations are, however, not cor- 
rected. Nevertheless, the characteristic features of the data 
are excellently reproduced. We see the two bump struc- 
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Fig. 10. Invariant differential cross section for the production 
of p,p,-K^ ,n~ , K'^ and K~ as a function of Feynman's vari- 
able Xf- The symbols denote the data at various values 
of the transverse momentum pt, histograms are ROC results. 
Data and calculations are multiplied by the factors given in 
the legend. 



ture at low and the shrinkage of the distributions with 
increasing multiplicities as well as their broadening with 
increasing energy. In the ROC calculations this behavior 
comes from the decay of the leading FBs situated near 
the projectile and the target rapidity at low multiplici- 
ties, while with increasing multiplicity the contributions 
from the increasing number of FBs dominate and due to 
momentum conservation the mean rapidity of the leading 
FBs change to lower values too. 

The same data were described by ref. with a sim- 
ple cluster model and by |Q in the framework of the 
dual parton model ||l[. Both attempts failed in reproduc- 
ing the two-bump structure at low multiplicities, obviously 
due to the absence of a diffractive component in these ap- 
proaches. 

In fig. ^ transverse momentum distributions for vari- 
ous particles are compared with experimental data. Both 
data and model results show roughly an exponential be- 
havior. The calculated spectra result from a convolution 
of the transverse momenta of the FBs with the internal 
momentum distributions of the primary particles in the 
FBs and the momentum distributions of secondaries in 
which primary resonances decay. At the highest momenta 
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an underestimation of the measured cross section is ob- 
served while the overall agreement is quite satisfactory 
for all particle types. 

Proton spectra as a function of a longitudinal variable 
like Feynman's variable x p are of special interest because 
of the peak from diffractive scattering for xj? — > 1. As 
already mentioned there is no special parameter which 
forces this channel to have a definite probability. There- 
fore, the good overall reproduction of the proton spec- 
tra in fig. H is quite remarkable, although the height of 
the peak is underestimated by the calculations. A better 
agreement might be obtained by giving the leading FBs a 
lower temperature than the central ones in analogy to the 
spectator-participant picture of hadron-nucleus reactions. 
This ansatz will be left to forthcoming considerations. 

Also at lower energies excellent agreement between 
data and ROC model calculations is achieved, as demon- 
strated in fig. 
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There the dependence of invariant dif- 
ferential cross sections on Feynman's variable xp for the 
production of various hadron species at a laboratory mo- 
mentum of 100 GeV/c is depicted. While the proton spec- 
tra show pronounced diffractive peaks for a;_F — > 1, the 
cross sections for the other particle types strongly decrease 



> 20 



if) 



X) 
C 
"D 
\ 

b 

X) 







10 




10 






90" 



135" 





60° 



105" 



150" 



75° 




120" 



0.2 



0.4 



0.2 0.40 0.2 0.4 

kinetic energy T(GeV) 

Fig. 12. Differential cross section as a function of the kinetic 
energy for vr^ production at laboratory angles between 15° and 
150°. Dots are data |37[, histograms ROC results. 



at large xp- The squared longitudinal matrix element 
with the parameter /3 in eq.(P) governs this behavior. It 
keeps the leading FBs at large \xp\ and the particles from 
the central FBs at small \xp\. The smaller (3 is chosen the 
larger the cross sections become for mesons and p's for 
Xp ^ 1. 

In stepping down the energy scale let us give a further 
example in fig. Shown are the invariant cross sections 
for three types of baryons as a function oi x — 2p;/y^ (pi 
being the longitudinal momentum) and for charged pions 
as a function of the rapidity . The main features of the data 
are well reproduced, although again the diffractive peaks 
in the proton spectra are underestimated. In case of the 
neutron spectra there are deviations especially at higher 
values of pt, while the A spectra are quite well described. 
Also the rapidity plateau at small and its shortening with 
increasing transverse momenta in the pion spectra is well 
reproduced. 

Finally, we verify the applicability of the ROC model 
at energies below the production threshold for strange par- 
ticles and show the kinetic energy spectra of 7r+ mesons 



at an incidence energy of 740 MeV in fig. 12. The peak at 
the high-energy end of the calculated spectra comes from 
the binary reaction pp — > dTr"*" . Due to insufficient energy 
resolution this peak is not seen in the measured spectra. 
The data are well reproduced in forward and backward 
direction while deviations become noticeable at sidewards 
angles. At low energy the parameter values discussed in 
sect. ^ are no longer important. Instead, the details of 
the treatment of resonance decays play an important role. 
Most of the 7r+ meson production proceeds via the cre- 
ation of A resonances. The spectra are therefore strongly 
influenced by the mass distributions of the decayed res- 
onances. In the present version of the ROC model the 
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current masses of all types of resonances are sampled by 
using a constant width of the Breit-Wigner distribution. 
This might be the reason for the observed deviations. 



3.3 Correlations 

The data considered so far can obviously be quite well re- 
produced by both string and FB models. In this subsection 
the question is discussed whether rapidity correlations are 
more sensitive with respect to their interpretation in terms 
of strings or FBs. The existence of short-range correlations 
is experimentally well established and many papers deal 
with the various aspects of correlations (see, e.g., ||38|,p9|, 
|^,|l|,|4|||). 

When looking for rapidity correlations one usually de- 
fines a single-particle 

Pi{y) = c^in^dcr/dy 
and a two-particle rapidity density 

P2{yi,y2) = (TiV/d^ff/dyidya- 

The latter is proportional to the probability of finding 
one particle at yi and a second one at 1/2 ■ In order to 
see whether the joint production of a pair of particles at 
(2/1 J 2/2) differs from an independent production of the two 
particles the two-particle correlation function 



p(200GeV/c)-Hp ^ p,(ly,l^0.25) + p2+X 



C{vi,v2) ^ P2{yi,y2) - pi{yi)pi{y2 



(19) 



is introduced as the difference between the two-particle 
density and the product of two single-particle densities. 
Non- vanishing values of C{yi, 1/2) indicate the presence of 
correlations. 

In fig. |l3| the measured and calculated two-particle 
correlation functions for different charge combinations of 
the two observed particles are compared. The pronounced 
peaks in the measured spectra are satisfactorily repro- 
duced by the ROC model calculations (histograms). How- 
ever, this fact alone is not yet a proof of the existence of 
FBs for at least two reasons. First, the measurements are 
inclusive ones and mixing of events with different mul- 
tiplicities can cause strong correlations as, e.g., pointed 
out in refs. Second, the presence of resonances 

among the emitted particles also tends to group the ob- 
served particles into clusters. That means, the observed 
correlation spectra contain always a superposition of ef- 
fects from resonances and other possible short-range phe- 
nomena. 

In order to see whether resonance production alone 
can reproduce the observed correlation patterns we carry 
out calculations with a special version of the ROC model 
(abbreviated by ROCS in the following). All FBs are com- 
pelled to degenerate into single hadrons by restricting the 
quark content of the FBs to be either qqq or qqq or qq. In 
this way all correlations implied by the presence of FBs 
are excluded. An additional parameter Wb is necessary in 
this case fixing the probability of baryon creation relative 




Fig. 13. The inclusive two-particle correlation function 
C(j/i,j/2) from 200GeV/c pp interactions as a function of 
i/2 — yi with fixed < 0.25 for different charge states of 
the two observed particles indicated by C'^'^, C and C' . 
Dots are data @, histograms ROC model, full lines PYTHIA 
and dashed lines ROCS results. 



to meson production, because the algorithm for the com- 
bined building of baryons and mesons from quarks (see 
sect. ^.2| ) is not applicable here. The value of W-q = 0.15 
adjusted to reproduce p production is similar to the prob- 
ability of diquark creation in striiig models, where, e. 5., 
the default values used in VENUS § and PYTHIA (§ are 
0.12 and 0.10, respectively. The ROCS version has been 
proven to reasonably reproduce most of the data discussed 
in the previous sections with readjusted parameters of 
[see eq. (^] and an increased radius parameter R. The 
ROCS results for the correlation are shown in fig. |l3| by 
the dashed lines, which underestimate the measured val- 
ues remarkably. We consider this result as a direct verifi- 
cation of the presence of short-range phenomena beyond 
resonance production. 

As a further proof of such additional short-range phe- 
nomena calculations with the string model PYTHIA Q| 
version 6.115 are carried out using the multiple interaction 
approach of ref. Q with varying impact parameter. The 
mean multiplicity of charged particles is reproduced by 
adjusting the regularization scale p_Lo (PARP(82)) of the 
transverse momentum spectrum. A Gaussian is taken for 
the matter distribution of the interacting protons, what 
gives a reasonable reproduction of the multiplicity distri- 
bution of charged particles at ^/s = 63 GeV. All other 
parameters are kept at their default values. The results, 
the full lines in fig. |l^, underestimate the data too. This 
finding confirms again that resonance production alone is 
insufficient for a correct reproduction of correlation data. 
In addition, the similarity of ROCS and PYTHIA curves 
shows that the fragmentation of multiple strings in PYTHIA 
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Fig. 14. Inclusive (n=all) and semi-inclusive two-particle cor- 
relations Cn(?7i,7?2) vs. 7)2 at fixed Tji — for the observed 
multiplicities n indicated in the figure from pp interactions at 
y/s — 23 GeV and 62 GeV. At 23 GeV the corresponding mean 
values of the true multiplicities are 9.1, 14.0 and 21.7, and the 
curves are calculated for multiplicities of 8-10, 14 and 18-24, 
respectively, while at 62 GeV the mean values are 8.1, 18.8 
and 30.8, and the curves are calculated for multiplicities of 8, 
18-20 and 28-34, respectively. Dots are data [^5|, histograms 
ROC model, full lines PYTHIA and dashed lines ROCS results. 



can be quite well imitated by creating hadrons in longitu- 
dinal phase-space. 

Correlations caused by the mixing of events with differ- 
ent multiplicities can be excluded by fixing the multiplic- 
ities of the considered events. This type of measurements 
will be called semi-inclusive in the following. For the cor- 
relation function the same formula (|l^) holds with the 
single and two-particle densities taken from events having 
a definite multiplicity. As an example results of Amen- 
dolia et al. are presented in fig. |l^ for the reaction 
p + p ci + C2 + X at two energies with ci , C2 and X 
standing for two charged particles and anything, respec- 
tively. The inclusive two-particle correlation (on top of 
the figure) is compared with semi-inclusive correlations at 
three narrow intervals of the observed multiplicities. The 
latter data sets, however, contain events whose true multi- 
plicities span a much wider interval (see fig. 5 in ref . ||45|| ) . 
Lacking the exact knowledge of the detector response we 
have selected multiplicity intervals around the mean val- 
ues of the true multiplicities in the calculations. Figure 
[1^ shows that the inclusive correlation function is much 
broader than the semi-inclusive ones, a feature of the data 
which is well reproduced by the ROC calculations. The 
semi-inclusive data with their distinctly smaller peaks and 
and pronounced dips of both sides of the peaks are well 
described too. Alone the peaks in the data for the high- 
est multiplicity especially at y/s — 23 GeV are underesti- 
mated. The PYTHIA and ROCS curves clearly underesti- 
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Fig. 15. The same as fig. y, but for semi-inclusive two- 
particle correlations C„(»7i, 772jvs. 772 at fixed ryi = —0.87. 



mate inclusive as well as semi-inclusive data. Again they 
are in good agreement with each other. 

In fig. |l5| the correlation function with one charged 
particle fixed outside midrapidity is given. The data are 
well reproduced by the ROC model, except for the high- 
est multiplicity interval at -^/s = 23 GeV. These deviations 
become larger if the distance of the fixed particle from 
midrapidity is increased (see fig. |l|). The ROCS version 
underestimates also in the non-symmetric cases all con- 
sidered correlations and agrees with the PYTHIA results 
with the exception of the lowest multiplicity interval at 
Vs = 62 GeV. There, the PYTHIA curve is closer to the 
experiment than the ROCS result, but still far from a good 
reproduction of the data. 

To better understand the results shown in figs. 
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the contributions from different subprocesses are de- 
picted in fig. ^ separately. In PYTHIA diffractive and 
nondiffractive processes are treated in different ways, while 
in the ROC model a diffractive process is simply one of 
the possible final channels without any special assump- 
tions. Consequently, the ROC results for all subprocesses 
exhibit the typical correlation pattern. On the other hand, 
the PYTHIA results show a large variety of completely 
different correlation functions in dependence on the con- 
sidered subprocess and on the pseudo-rapidity interval 
where one of the particles is fixed. In the first row of 
fig. [1^ the diffractive excitation of the projectile is con- 
sidered. We find the target proton after the interaction at 
negative rapidities and the excited projectile as a string 
in case of PYTHIA and as a number of FBs in case of 
ROC calculations mainly at positive rapidities. Values of 
the correlation function around zero near the target ra- 
pidity indicate that the emission of the target proton is 
weakly correlated with the emission of the charged par- 
ticles from the excited subsystem(s) (see the symmetric 
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Fig. 17. Semi-inclusive two-particle correlations C„ (771, 772) 
vs. 772 at several fixed values of 771 indicated in the figure for 
the multiplicity interval = 10 — 11 at ^/s — 62GeV. The 
data (dots) are compared with ROC model (histograms) 
and PYTHIA (full lines) results for single- {pp — > Xp and 
pp — > pX), double- {pp pXp) and nondiffractive {pp XX) 
processes. 



case 771 =0). For the diffractive excitation of the target 
(the second row) the unchanged target proton is at pos- 
itive, the excited subsystems at negative rapidities. The 
whole picture is simply reversed in the symmetric case 
(771 = 0), while we get from PYTHIA completely differ- 
ent pictures for the non-symmetric measurements with 
771 = -0.87 and 771 = -1.96. The PYTHIA result for 
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Fig. 18. Number of charged particles rich versus number of 
FBs A'^ in ROC model calculations for p+p at = 23 GeV and 
62 GeV. The size of the boxes is proportional to the calculated 
cross section. The full line in the upper part represents the 
mean value of charged particles (rich) and the line in the lower 
part the FB size (rich) /A as a function of N. 



double diffraction (the third row) is of special interest. 
Both interaction partners remain unchanged and the ex- 
cited subsystem is mainly at central rapidities. The central 
string acts like a cluster or FB and the results from ROC 
and PYTHIA arc similar. The contribution of this special 
process is responsible for the somewhat better reproduc- 
tion of the data by PYTHIA in the multiplicity interval 
n = 10—11 at = 62 GeV in fig. |l^. With increasing 
multiplicities the diffractive contributions become smaller 
and the result is dominated by the nondiffractive compo- 
nent. In this case (the fourth row) we have in PYTHIA 
the typical multi-string picture with a clear underestima- 
tion of the measured correlation. Altogether we see that 
even at fixed multiplicities the contributions from differ- 
ent processes may result in distinct diffraction patterns in 
a model dependent way. 

In spite of the results presented here it should be men- 
tioned that in the literature examples for the description 
of short-range correlations bystring models can be found. 
So the dual parton model |l|, which has a multi-string 
structure similar to PYTHIA, reproduces correlation data 
either under certain assumptions for the fragmentation 
functions (see, e.g., [^,|4^) or as the result of the decay 
of minijets and of soft strings with sea quarks at their 
ends Ipn]. Such strings emit hadrons into limited rapidity 
regions like clusters. 

In order to clarify the origin of the deviations be- 
tween ROC calculations and data at high multiplicities 
(see fig. |l6|) we show the correlation between the num- 
ber of FBs and the number of charged particles in fig. ^ 
First of all we observe large fiuctuations regarding the FB 
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size, where FB size is understood as the mean number of 
charged particles originating in a FB. At ^/s = 62 GeV 
it is possible that, e.g., 18 charged particles may be pro- 
duced in events with numbers of FBs ranging from 2 up 
to 20. The mean number of charged particles (the line 
in the upper part of fig. Iq) flattens out for large num- 
bers of FBs especially at -/s = 23 GeV. That means that, 
although we use constant parameters and R, the FB 



Vs=24GeV 45GeV 



63GeV 



size (the histogram in the lower part of fig. |18) becomes 
smaller with increasing number of FBs due to the phase- 
space factor, and this influence is at ^/s = 23 GeV much 
stronger than at 62 GeV. Here we see the origin for the 
discrepancy with the data at high multiplicities. This is 
confirmed by a series of calculations with increased pa- 



rameter values of R (and accordingly adapted Af , eq. (| 
to keep the mean number of charged particles constant). 
Due to the larger R the size of the FBs is increased and 
a good reproduction of the data at high multiplicities is 
possible while the low-multiplicity data are overestimated. 
That means, the present version of the ROC model is too 
simple with regard to the FB size. The volume V derived 
from the parameter R is considered as a measure for the 
overlap region of the two colliding hadrons which defines 
the spatial region where particle production takes place. 
Therefore, it is understandable that a value of V inde- 
pendent of the impact parameter is surely an oversimpli- 
fication. Since the number of produce d F Bs is correlated 
with the impact parameter (see sect. 2.1), an increase of 
R with increasing number of FBs would be quite a nat- 
ural improvement of the model. A dependence of the FB 
temperature and of the momentum distribution of FBs 
described by P and 7 on the number of FBs cannot be 
excluded as well. A central collision may produce hotter 
FBs with a more isotropic phase-space distribution than 
a peripheral one. Such a view is supported by an analysis 
[ pl| of long-range correlations between charged particles 
emitted into the forward and backward hemispheres at 
much higher energies (Vs 0.3 - 1.8 TeV). The E735 col- 
laboration came to the conclusion, that the cluster size 
may increase as a function not only of ^/s but also of the 
particle multiplicity. 

An analysis of long-rang correlations in the energy re- 

fa considered here has been carried out by Uhlig et al. 
. Their results are shown in fig. |l^. The mean number 
of charged particles emitted into the backward hemisphere 
is plotted as a function of the number of forward emitted 
particles. The slope of these curves is a measure for the 
strength of the correlation between particles ejected into 
the various regions of the hemispheres. Although there 
are some systematic deviations at low multiplicities and at 
yfs — 24 GeV the overall trend of the data is well described 
by the calculations. We see the strongest correlation for 
adjacent regions in pseudo-rapidity 77 while the strength 
is smaller if there is a gap between the considered 77 inter- 
vals. The emission is nearly independent for the regions 
in the lowest row of fig. ^ which are outside the central 
rapidity plateau. The authors of ref. explain their re- 
sults with the correlation between clusters consisting of 
about three particles (neutral ones included). In fig. ^9] 
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Fig. 19. Mean number of charged particles emitted into the 
backward hemisphere (ns) as a function of the number of 
charged particles emitted into the forward hemisphere np for 
symmetric intervals of the pseudo-rapidity r] indicated in the 
figure, where /(s) = 0.5 • ln(s ■ GeV"^) - 2. Dots are data H, 
histograms ROC model, full lines PYTHIA and dashed lines 
ROCS results. 



we see a striking agreement between the results of ROC, 
ROCS and PHYTHIA calculations. Obviously, both FB as 
well as string models are able to reproduce the observed 
long-range correlations in the considered energy region. 



4 Conclusions 

We have presented the empirical ROC model for soft hadron 
production. It is based on the parton picture of hadrons as 
well as on statistical and thermodynamical considerations. 
Experimental results from pp interactions in the energy re- 
gion between particle production thresholds and ISR en- 
ergies can be well described with a moderate number of 
parameters which are either constant or have a smooth en- 
ergy dependence. A comparison of the ROC fireball model 
and the PYTHIA string model with regard to short-range 
rapidity correlations seems to favor the ROC model, al- 
though the description of short-range correlations is pos- 
sible in string models too. Thus we do not have a clear 
answer to the problem of hadron production via strings or 
fireballs. A systematic consideration of all available data 
especially those at higher energies is necessary to come to 
a decision between the possible scenarios: strings, fireballs 
or something in between. 
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